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Abstract 

Near the critical point, isothermal interfacial zones are investigated 
starting from a non-local density of energy. From the equations of mo- 
tion of thermocapillary fluids pQ , we point out a new kind of adiabatic 
waves propagating along the interfacial layers. The waves are asso- 
ciated with the second derivatives of densities and propagate with a 
celerity depending on the proximity of the critical point. 

Resume 

Pres du point critique, les couches interfaciales sont modelisees a Vaide 
d 'une densite d 'energie non locale. A partir des equations du mouvement des 
fluides thermocapillaires [lj, nous mettons en evidence des ondes adiabatiques 
se propageant le long des couches interfaciales. Ces ondes associees aux 
derivees secondes des densites se meuvent avec une celerite dependant de la 
proximite du point critique. 



Version francaise abregee 

Le modele le plus simple permettant de considerer les couches capillaires 
et les phases comme un unique milieu continu, consiste a prendre en compte 
une energie comme la somme de deux termes : le premier correspond a 
l'energie du milieu suppose uniforme et de composition egale a la composi- 
tion locale, le second est associe a la non uniformite du fluide et exprime par 
un developpement en gradient de la masse volumique qui est limite au second 
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ordre [21 E]- Ce modele permet le prolongement au cas dynamique des etudes 
effectuees sur les systemes en equilibre. L'energie interne volumique du flu- 
ide est maintenant proposee sous la forme d'une densite a dependant non 
seulement de grad p mais aussi de grad s ou p et s notent respectivement la 
masse volumique et l'entropie specifique. Les milieux associes appeles fluides 
thermocapillaires [JJ ont une densite d'energie de la forme (pQ). Pour un flu- 
ide isotrope, l'energie prend la forme ([2]) et l'equation du mouvement s'ecrit 
sous les formes equivalentes ((SJ) ou (jSJ) dans lesquelles une entropie h et une 
temperature 9 dites thermocapillaires explicitees en (j3J) font intervenir deux 
nouveaux vecteurs <fr and ^ associes a la non homogeneite des zones interfa- 
ciales et donnes par Une telle energie interne associe dans les equations 
du mouvement et de l'energie un tenseur des contraintes spherique. Le cas 
des ecoulements isentropiques associes a l'equation (J7J) correspond a la con- 
servation de l'energie (jSJ). Les mouvements isothermes represented par (Q et 
fTTOj) ont comme cas particulier les equilibres isothermes entre phases. 
II est alors possible d'etudier les profils des densites dans les equilibres uni- 
dimensionnels correspondant aux zones interfaciales planes. Le choix d'une 
energie interne thermocapillaire de la forme (|lip associee a une energie in- 
terne a du milieu suppose homogene ramene l'etude des zones interfaciales 
a l'analyse d'un systeme dynamique (112bj) . (112aj) qui donne deux equations 
differentielles. L'analyse asymptotique de ce systeme au voisinage du point 
critique montre que la densite de masse pilote le systeme en ce sens que Ton 
est ramene a une decomposition dynamique representee par le systeme (115bj) . 
(j 15 aft . L'integration de ce systeme nous ramene simplement au cas d'un 
modele ne faisant plus intervenir que grad p. On retrouve immediatement 
dans (fTEj) . les resultats proposes par Rowlinson et Widom [7] pour le profil 
de la masse volumique, l'epaisseur de la couche interfaciale et la valeur de la 
tension superficielle des interfaces fluides au voisinage du point critique. 
II est maintenant possible d'etudier, dans les interfaces isothermes, les on- 
des d'acceleration correspondant a des discontinuites faibles des mouvements 
isentropiques c'est-a-dire celles pour lequelles p,s, gradp et grads sont con- 
tinus a la traversee des surfaces d'onde. Un calcul a la Hadamard [TT] est 
propose. II prend en compte l'equation de conservation de la masse (ITS]) , 
l'equation du mouvement sous la forme (jSJ) et la condition de Rankine- 
Hugoniot associee a l'equation du mouvement ecrite sous la forme equivalente 
(jSJ). II permet alors d'obtenir un systeme fl23l) . fT25|) . fl27|) de trois equations 
lineaires et homogenes relatif aux discontinuites des derivees secondes a 
travers les surfaces d'onde (paragraphe 4.1). La condition de compatibilite 
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de ce systeme s'exprime par l'equation (I2"gj) qui donne ainsi la celerite de 
l'onde. 

La relation (I15al) issue du systeme dynamique fll4bj) . (I14al) obtenu par l'ana- 
lyse asymptotique du paragraphe 3, permet d'expliciter la relation liant les 
gradients d'entropie et de masse volumique. Nous la calculons pour la valeur 
de la masse volumique critique. Nous en deduisons la valeur explicite de la 
celerite des ondes d'acceleration cisaillees dans les interfaces fluides. Cette 
celerite exprimee par la relation (I29I) depend des conditions thermodyna- 
miques du point critique et de l'importance des coefficients associes aux ter- 
mes representant l'inhomogeneite du fluide dans la zone interfaciale. Elle est 
proportionnelle a (T C — T Q ) 2 ou T represente la valeur de la temperature dans 
les phases liquide et vapeur et T c note la temperature critique. 
Les ondes solitaires dans la direction normale aux zones interfaciales ne 
dependent pas du gradient d'entropie. II n'en est pas de meme des ondes 
isentropiques d'acceleration qui se meuvent le long des interfaces : le fait que 
l'energie interne depend non seulement du gradient de masse volumique mais 
aussi du gradient d'entropie fait apparaitre un nouveau type d'onde qu'il est 
impossible de mettre en evidence dans les modeles plus simples ne faisant 
intervenir que le gradient de masse volumique. Ces ondes exceptionnelles au 
sens de Lax [12] necessitent le cadre d'un modele physique a au moins deux 
dimensions d'espace. Des experiences effectuees recemment dans des labora- 
toires embarques dans des engins spatiaux et utilisant comme fluide test le 
gaz carbonique dans des conditions qui le placent au voisinage de son point 
critique semblent montrer l'existence de telles ondes [13] et justifieraient alors 
le bien fonde du modele de fluides thermocapillaires. 

1 Introduction 

To study capillary layers and bulk phases, the simplest model considers an 
energy as the sum of two terms: a first one corresponding to a medium with 
a uniform composition equal to the local one and a second one associated 
with the non- uniformity of the fluid [21 E]- The second term is approximated 
by a gradient expansion, typically truncated to the second order. A repre- 
sentation of the energy near the critical point therefore allows the study of 
interfaces of non-molecular size. Obviously, the model is simpler than mod- 
els associated with the renormalization-group theory [4] . Nevertheless, it has 
the advantage of extending easily well-known results for equilibrium cases 
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to the dynamics of interfaces [3,E]- For equilibrium, Rowlinson and Widom 
[7] pointed out that the model can be extended by taking into account not 
only the strong variations of matter density through the interfacial layer but 
also the strong variations of entropy In dynamics, for an extended Cahn 
and Hilliard fluid, the volumic internal energy e is proposed with a gradient 
expansion depending not only on grad p but also on grad s (p is the matter 
density and s the specific entropy) : 

£ = f(p, s,grad p,grad s). (1) 

The medium is then called a thermocapillary fluid [T] . Using an energy in 
the form flTJ, we have obtained the equations of conservative motions for 
nonhomogeneous fluid near its critical point Pfl El E] • 

The idea of studying interface motions as localized travelling waves in a 
multi-gradient theory is not new and can be traced throughout many prob- 
lems of condensed matter and phase-transition physics [10J. In Cahn and 
Hilliard's model [3J, the direction of solitary waves was along the gradient 
of density [5J [TU]. Here, adiabatic waves are considered and a new kind of 
waves appears. The waves are associated with the spatial second derivatives 
of entropy and matter density. For this new kind of adiabatic waves, the 
direction of propagation is normal to the gradient of densities. In the case of 
a thick interface, the waves are tangential to the interface and the wave celer- 
ity is expressed depending on thermodynamic conditions at the critical point. 



2 Equations of thermocapillary fluid motions 

The equations of motion are proposed in [H [HI IS] - Due to the fact the fluid is 
isotropic, grad p and grad s are taken into account by their scalar products 
only. Let us denote 

(3 = (grad p) 2 , x — grad p . grad s, 7 = (grad s) 2 . 

In variables p, s, (3, x, 7, 

£ = 9(p,s,/3,x,l)- (2) 

The equations of thermocapillary fluids introduced two new vectors $ and 
^ such that : 

<fr = C grad p + D grad s, & — D grad p + E grad s, (3) 
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with 

C = 2 £ p : D = £ jX , E = 2 e y. 

We denote 

h = e )P - div <I>, 6 1 = - (e s - div . (4) 

In the case of compressible fluids, scalars e p and (l/p)e s are the specific 
enthalpy and the Kelvin temperature. Look at two particular cases : 

2.1 Conservative motions 

We obtained [I] the equation of motion in the form : 

r = 9 grads - grad(/i + Cl), (5) 

where T is the acceleration vector, Q is the extraneous force potential. This 
equation is equivalent to the balance of momentum: 

d 

— (pu) + div(pu ® u — a) + pgradfi =0 (6) 

with <7^ = — (P — p div $) 51 — & p { — \E , - J S j where P — p£ tP — £ (in the case 
of classical compressible fluids, P denotes the pressure) and u is the fluid 
velocity. For conservative motions, 

which is equivalent to the balance of energy 

de dQ 

— +div((e-(r)u) -divW-p— = 0, (8) 

with W = $ + ^ * and e = | pu 2 + e + pVt. 

at at z 

2.2 Isothermal motions 

We obtained [1] the equation of motion in the form : 

e = T , (9) 

r = - grad(p + fi), (10) 

where T is constant and p = £ )P — sT — div $ is the chemical potential of 
the thermocapillary fluid. 
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3 Liquid- vapor interface near its critical point 



The critical point associated with the equilibrium of two bulks of a fluid 
corresponds to the limit of their coexistence. The thickness of the interface 
increases as the critical point is approached and it becomes infinite when the 
interface itself disappears. As its critical point is approached, the gradients 
of densities are then smooth. In the following, we consider the case when 

e = pa(p, s) + - ((J (grad p) 2 + 2D grad p . grad s + E (grad s) 2 j (11) 

where a denotes the specific internal energy of the fluid in uniform com- 
position (let a(p, s) be the analytic a as it might be given by a mean-field 
theory), C, D, E are constants and CE — D 2 > 0. If D = and E = 0, we 
are back to the Cahn and Hilliard model of capillarity [3]. If not, we deduce 

h = h -(C Ap + D As), 6 = T- - (DA p + E As), 

P 

where ho = a + p a'„ and T = a' s are respectively the specific enthalpy and 
the Kelvin temperature of the homogeneous fluid of matter density p and 
specific entropy s; we denote by p = e' p — s T its chemical potential. 
At phase equilibrium, Eq. ([9]) is verified when To is the temperature in the 
bulks. If we neglect the body forces, we obtain 

DAp + EAs = e' s -pT , (12a) 
CAp + DAs = e' p -sT -p 1 (12b) 

where p l is the value of p in the liquid and the vapor bulks. 

To be in accordance with Rowlinson and Widom ([7J p. 253), we consider a 

representation of pct(p, s) in the form 

B ( / \ 2 / A 2 \ 2 \ 

pa(p,s) = — i[A(p- p c ) 2 + psj + [ps + — (T c - T ) J J + pp l + psT 

(13) 

where A and B are two positive constants associated with the critical condi- 
tions and s which is defined except to any additive constant is choosen such 
that s c = 0. System (12) yields 

D Ap + E As = ^p(p-p c ) 2 + 2^-p 2 s + p(T c -T ), 
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C Ap + D As = 2B (p - p c ) 3 + 2^ps (p - p c ) + ^(p - p c ) 2 s 
+ 2^ps 2 + s(T c -T ). 

The equations of the coexistence curve at T = To are associated with A p = 

A 2 

andAs = 0. They are as in [7] , A(p— p c ) 2 +ps = and ps + — (T c —T ) = 0. 

B 

In the bulks, near the critical point, the respective magnitudes of p — p c and 
s with respect to T c — T are p — p c ~ cte (T c — T ) 5 , s ~ cte (T c — T ). In 
one-dimensional problems, p = p(y), s = s(y). To consider the respective 
order of magnitude of the densities and the physical scales associated with 
the interfacial sizes, we look at the change of variables 

T c -T = eT, y = e n Y, p(y) - p c = R(Y), s{y)=eS(Y), 

where 0<e<l, and n is constant. The main part of system (12) leads to 

2 „ / i d 2 R „ d 2 S\ ( B 2 B - „ \ . A . 

-2n ( d 2 R d 2 S \ 3 f g „ B 



Ct ^ + Dt W^) = ^Bi^-pflS) (14b) 

Then n — —\ and in this asymptotic analysis, if " denotes the second deriva- 
tive with respect to the space variable y, system (12) yields the approximation 

2ps = ^-(T -T c )-A(p-p c ) 2 , (15a) 
Cp" = 2B(p-p c ) 3 + 2jps(p- Pc ) (15b) 

and consequently by elimination of s between ( 115bl) and (1 15 aft , we obtain 

C p" = B (p - Pc f -A(T C - T )(p - p c ), (16) 
and for planar liquid-vapor interfaces, 

1 




2 C P' 2= I ^l>-l>r) 2 -^h«r--J\>) ) ^ (17) 

We are back to the classical van der Waals theory for mass density profile in 
the interfacial layer (see [7J, p.p. 250-251). 
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4 Weak discontinuity in conservative motions 



4.1 Conditions of a weak discontinuity 



In a fixed coordinate system, to describe the fluid motion, we refer to the 
coordinates x — ( 00 ■ 00 ■ 00 ~\ cLS the particle position (Eulerian variables). We 

the time-space variables; V 



denote by z 



u 



is the time-space 



N 



N*V = n*u — g is the fluid velocity with respect to S t 



velocity. Let us consider a mobile surface S t defined in the physical space 

occupied by the fluid, we denote by g the celerity of E 4 , n its normal vector, 

-9 
n 

(—v is the celerity of E t with respect to the fluid) and n* is the transposed 
of n. 

To Eq. (Q of motion, we have to add the equation of balance of mass 

dp 



at 



+ div(pu) = 0, 



(18) 



Eq. ([7]) of conservation of the entropy and relation ([T3|) . 
Due to the important fact that the first discontinuities of p and s are not ex- 
isting in thermocapillary mixtures, weak discontinuities of isentropic motions 
correspond to p, s, dp/dx, ds/dx, continuous through the wave surfaces. As 
in [11], we denote by [ ] the jump of a tensorial quantity through a surface of 
discontinuity S t . The jump of a spatial derivative of a continuous quantity is 
colinear to the normal vector to the wave surface; as in Hadamard's tensorial 
framework, there exits two scalar fields Ax, A2 on S t , such that 



dp 
dz 
ds 
dz 

and consequently 











9_,dfK* 
dz K dz> 
d ds * 



AiNN* 
: A 2 NN* 



d_ ,dp_. 
<9x W 
d ,ds s 
dx"fe 



From [V] 


H 



Xi-- 

0, we deduce 



[Ap] and \ 2 = [As~ 
dV 



dz 



E N* with N* 



: Ainn*, (19a) 
= A 2 nn*, (19b) 

(20) 

-g, n*) and S = 



where H is a 3- vector field on £+. Then, 



du 

dz 



V = HN* V = f H. 
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Equation of mass conservation (TTgj) is equivalent to 

dp 



dz 



V + pTr 



dV 
dz 



where Tr denotes the trace operator. Then, 
quently, 

n* H = 0. 

Eq. (JTj) of conservation of entropy implies 



0. 



Tr 



dV\ 
~dz~) 



and conse- 
(21) 



d ds, 
dz~^di' 







d ( ds 



dz \dz 



0. 



Then, 



dz \dz 



dz J \ dz J 



d / ds \ * 

Due to the fact that — ( — ) is a symmetric tensor, 

dz\dzJ 



\dz 



A 2 NN*V + N3* — =0 <£> N [ \ 2 v + —H 



0. 



ds 



<9x 



or, 



ds 



X 2 v + — H = 0. 
ax 



(22) 



From Rankine-Hugoniot condition associated to Eq. ([6]), we obtain the com- 
patibility condition: [— gpu* + pn*uu* — n*a] = 0, and the continuity of 
p, dp/dx, s, ds/dx, yields [div <&] = or [CAp + DAs] = which is equiv- 
alent to 

CAi + DX 2 = 0. 



Consequently, there exists a scalar field A3 on Ej such that 

d 



(CAp + DAs) 



A,n*. 



. <9x 

Equation of motion ([5]) yields 

[r] = [#]grad S -A 3 n, 
or, p v H = — (DXi + E\ 2 ) grad s — p A3 n. 



(23) 



(24) 



9 



By projection on the normal and tangent plane to T, t and taking relation 
( 12 ip into account, we obtain 

D (n* grads) \ x + E (n* grads) A 2 + p A 3 = 0, , (25) 

p v H = - (D Ai + E A 2 ) grades, (26) 

where grad tg s denotes the tangential part of grads in £ t . Elimination of H 
in the relation (126]) comes from relation (1221) . and we get 

D (grad tg s) 2 Ai + (E (grad tg s) 2 - pv 2 ) A 2 = 0. (27) 

Consequently, we obtain the system (1231) . (1251) . (127]) of three homogeneous 
linear equations with respect to the variables Ai, A 2 , A3. The compatibility of 
the three equations yields 

^ = (gg-g'Hg«W (28) 

4.2 Celerity of isentropic waves of acceleration 

The temperature in liquid and vapor bulks is To. Then relation (I15al) yields 

grad ( A(p — p c ) 2 + 2 p s ) = ^ A(p — p c ) grad p + grad ps = 0. 

The value p = p c of the matter density in the interface corresponds to the 
maximum value of grad p. The matter density p c is characteristic of the 
interfacial matter. For such a value, grad ps = and s = (A 2 /2Bp c ) (Tq—T c ). 

A 2 

2Bpl 

2 _{CE- D 2 ) A 4 (grad tg p) 2 (T c - T ) 



Consequently, grad s = - — - (T c — T ) grad p, and 



ACB 2 p 



A 2 

Due to relation ( 1TT1) . we obtain when p = p c , C (grad tg p) 2 = — — (T c — T ) 2 , 

2 D 

and consequently, 

2 (CE — D 2 ) A 6 (T c — T ) 4 . . 

" = ~ s6*p& ■ (29) 

In the interfacial layer, grad p is normal to iso-density surfaces. The isen- 
tropic waves of acceleration associated with a weak discontinuity shear the 
interfacial layer. The wave celerity which is proportional to (T c — T ) 2 , van- 
ishes at the critical point and can be calculated numerically by means of a 
state equation. 
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5 Results and discussion 



Near the critical point, for thermocapillary fluids, the variation of matter 
density leads to the variations of entropy (see Eq. (j 15 aft ). This extension of 
Cahn and Hilliard fluids seems at first sight purely formal and at equilibrium 
yields same results as the classical van der Waals model does [2] . In dynamics, 
motions in the normal direction to fluid interfaces (as solitary waves) are not 
involved in an additive dependance of the entropy gradient. In this paper, 
we see that the dependance of entropy gradient is necessary for isentropic 
waves of acceleration along the interfaces: the fact that the internal energy 
depends not only on the gradient of matter density but also on the gradient 
of entropy, yields a new kind of waves which does not appear in the simpler 
models. It is easy to see they are exceptional waves in the sense of Lax [12] 
and they appear only in, at least, two-dimension spaces. Recent experiments 
in space laboratories, for carbonic dioxide near its critical point have showed 
the possibility of such waves [T3] and should justify the well-founded of the 
thermocapillary fluid model. 
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